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Abstract 

1 /2 3/2 

Subleading Isgur-Wise form factor t(v ■ v ) at 0(l/m Q ) for A b — > A c [ ' ' 
weak transition is calculated by using the QCD sum rules in the framework 
of the heavy quark effective theory (HQET), where A c ( and A c [ are the or- 
bitally excited charmed baryon doublet with J p = (l~/2, 3 _ /2). We consider 
the subleading contributions from the weak current matching in the HQET. 
The interpolating currents with transverse covariant derivative are adopted 
for A^ 2 and A^( 2 in the analysis. The slope parameter p 2 in linear approxima- 
tion of r is obtained to be p 2 = 2.76 and the interception to be r(l) = —1.27 
GeV. 
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I. INTRODUCTION 



The ground state bottom baryon A b weak decays [1] provide a testing ground for the 
standard model (SM). They reveal some important features of the physics of bottom quark. 
The experimental data on these decays have been accumulated to wait for reliable theoretical 
calculations. With the discovery of the orbitally excited charmed baryons A c (2593) and 
A c (2625) [2], it would be of great interest for one to investigate the A b semileptonic decays 
into these baryons. 

From the phenomenological point of view, these semileptonic transitions are interesting 
since in principle they may account for a sizeable fraction of the inclusive semileptonic rate 
of A b decay. In addition, the properties of excited baryons have attracted attention in recent 
years. Investigation on them will extend our ability in the application of QCD. It can also 
help us foresee any other excited heavy baryons that have not been discovered yet. 

The heavy quark symmetry [3] is a useful tool to classify the hadronic spectroscopy 
containing a heavy quark Q. In the infinite mass limit, the spin and parity of the heavy 
quark and that of the light degrees of freedom are separately conserved. Coupling the spin 
of light degrees of freedom jg with the spin of heavy quark sq = 1/2 yields a doublet with 
total spin J = j e ± 1/2 (or a singlet if jg = 0). This classification can be applied to the 
Ag-type baryons. For the charmed baryons the ground state A c contains light degrees of 
freedom with spin-parity j[ = + , being a singlet. The excited states with jf = 1~ are spin 
symmetry doublet with J p (l"/2,3~/2). The lowest states of such excited charmed states, 
A*{ 2 and A 3 C { 2 , have been observed to be identified with A c (2593) and A c (2625) respectively 
[2]- 

However, the difficulties in the SM calculations are mainly due to the poor understanding 
of the nonperturbative aspects of the strong interaction (QCD). The heavy quark effective 
theory (HQET) based on the heavy quark symmetry provides a model-independent method 
for analyzing heavy hadrons containing a single heavy quark [3]. It allows us to expand 
the physical quantity in powers of 1/mg systematically, where uiq is the heavy quark mass. 
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Within this framework, the classification of the A b exclusive weak decay form factors has 
been greatly simplified. The decays such as A b — > A C W [4], A h — > Y^lv [5], A b — > e£"VZP 
[6], At — > p(A) [7] have been studied. 

To obtain detailed predictions for the hadrons, at this point, some nonperturbative QCD 
methods are also required. We have adopted QCD sum rules [8] in this work. QCD sum rule 
is a powerful nonperturbative method based on QCD. It takes into account the nontrivial 
QCD vacuum which is parametrized by various vacuum condensates in order to describe 
the nonperturbative nature. In QCD sum rule, hadronic observables can be calculated by 
evaluating two- or three-point correlation functions. The hadronic currents for constructing 
the correlation functions are expressed by the interpolating fields. In describing the excited 
heavy baryons, transverse covariant derivative is included in the interpolating field. The 
static properties of A b and A cl (A cl denotes the generic j[ = 1~ charmed state) have been 
studied with QCD sum rules in the HQET in Ref. [9] and Ref. [10,11], respectively. Recently, 
the leading order Isgur-Wise (IW) function is also calculated in the HQET QCD sum rule 
in Ref. [12]. 

In A b — > A cl decay, 1/mg corrections are very important. At the heavy quark limit of 
rriQ — > oo, the transition matrix elements should vanish at zero recoil since the light degrees 
of freedom change their configurations. Nonvanishing contribution to, say, B(A b — > A c \lv) 
at zero recoil appears at 1/ mq order. Since both A b and A cl are heavy enough, the behavior 
of the matrix elements near the zero recoil is very important. That explains why people 
pay attention to the next-to-leading order (NLO) contributions. The same situation occurs 
in heavy mesons. As for B — > D\{D*£)lv decay, leading and subleading Isgur-Wise (IW) 
functions have been computed using QCD sum rule in Ref. [13-17]. They showed that the 
branching ratio is enhanced considerably when the subleading contributions are included. 

In HQET, 1/toq corrections appear in a two-fold way. At the Lagrangian level, sub- 
leading terms are summarized in Ai and A2. Ai parametrizes the kinetic term of higher 
derivative, while A2 represents the chromomagnetic interaction which explicitly breaks the 
heavy quark spin symmetry. At the current level, 1/rriQ corrections come from the small 
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portion of the heavy quark fields which correspond to the virtual motion of the heavy quark. 
In this work, the subleading IW function from the latter case, i.e., at the current level, is 
analyzed in the HQET QCD sum rules. 

In Sec. II, the weak transition matrix elements are parametrized by the leading and 
subleading IW functions. By evaluating the three-point correlation function, we give the 
subleading IW function in Sec. III. We present, in Sec. IV, the numerical analysis and 
discussions. The summary is given in Sec. V. 

II. WEAK TRANSITION MATRIX ELEMENTS AND THE SUBLEADING 

ISGUR-WISE FUNCTIONS 



The weak transition matrix elements for A;, — > A cl are parametrized by the 14-form 
factors as 



m a c i(V, s ') f i7m + + -y 5 u Ab (v,s) , 
u\ cl (v', s') Gi7 M + G 2 v fM + G 3 v'^ u Ab (v, s) , 



(la) 
(lb) 



>/4M Acl(1/2) M Afc 

y/AM Acl{1/2) M Ab 

3 

{ -^7==^ - = *U*> s ') + K 2 v, + K^) + K 49 Ju Ab (v, s) , (1c) 

v /4JM Acl (3/2)MA b L J 

3 

(kl^v' s')\A,\k b {v^ = ^ g/) r {Ni ^ + N2V/i + NW ) + N4g J i5UAb{v ^ s) (ld) 

V 4M A c i(3/2)^A 6 L J 

where v (v') and s(s') are the four-velocity and spin of A b (A cl ), respectively. And the form 
factors Fi, G iy Ki and JVj are functions of y = v ■ v' . In the limit of mq — > oo, all the 
form factors are related to one independent universal form factor £(y) called Isgur-Wise 
(IW) function. A convenient way to evaluate hadronic matrix elements is by introducing 
interpolating fields in HQET developed in Ref. [18] to parametrize the matrix elements in 
Eqs. (1). With the aid of this method the matrix element can be written as [19] 



cTb = h^Th® = £(yKC,rv„ 



(2) 
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at leading order in 1/rriQ and a s , where T is any collection of 7-matrices. The ground state 
field, ip v , destroys the A b baryon with four- velocity v; the spinor field is given by 

C = € /2a + ^= (7° + v a )l^l /2 , (3) 

where ipl^ 2 is the ordinary Dirac spinor and ipv^ 2a is the spin 3/2 Rarita-Schwinger spinor, 
they destroy A^ 2 and Aj!{ 2 baryons with four-velocity v, respectively. To be explicit, 

Fi = -j={y - 1) Z(y) , G 1 = i=(y + l)^), 
f 2 = g 2 = ~ av) , k 1 = n 1 = t(y) , 

(others) = . (4) 

In general, the IW form factor is a decreasing function of the four velocity transfer y. Since 
the kinematically allowed region of y for heavy to heavy transition is very narrow around 
unity, 

1< y < A » Acl ~ 1.3 , 5 
- y - 2M Ai M Acl ' 1 J 

and hence it is convenient to approximate the IW function linearly as 

£(y) = £(i)(i-P 2 (y-i)), (6) 

where p| is the slope parameter which characterizes the shape of the leading IW function. 

The Aqcd/^q corrections come in two ways. One is from the subleading Lagrangian of 
the HQET while the other comes from the small portion of the heavy quark field to modify 
the effective currents. We only consider the latter case here. 

Including Aq CD /m 6 and Aqcd/^c, the weak current is given by 

Keeping the Lorentz structure, the subleading terms are expanded in general as 
htfifrh® = rAr[ c) v a i + r 2 (c W+ r 3 (c) 7.)rA„ , 

h$r0W = Wr(T®v a 1 + r 2 ( V/+ r 3 (6) 7 a)A, , (8) 



where are the subleading IW functions to be evaluated. 
The matrix elements of these currents modify Eq. (4) as 



= (y- 1)£ - e c \(y - 1)(- Tl (c) + r 2 (c) ) + 3r 3 (c) 



(y-lKr^-r^-r, 



-(b) 



V3F 2 
V3F 3 



= -2C + e c 



-2e, 



2y Tl (c) + 2r 2 (c) 



+ e& 



2rf ) + 2r 2 (6) 



(1 + y)r 2 w + r. 



V3Gi = (y + l)e - e c \(y + l)(r< c) + r 2 (c) ) + 3r 3 (c) 



+ e b 



-2£ + e c 



2yr[ c) + 2r 2 (c) 



2e, 



V3G 3 


= 2e b 


"(y 


-l)rf + r«; 








= £ + 




>) - r (c) " 

'1 '2 


+ 66 








= £- 


e c 


r (c) 4- r (c) " 

'1 + r 2 


+ £6 




+ r<«; 


K 2 


= iV 2 




-2e r (c) 









^3 
^4 



-N 3 = 2e 6 r 2 (b) , 



-iV 4 = 2e fe r 3 v 



(6) 



(9) 



where €q = l/2m<g. It is quite convenient to define 



^ r) = (7 Q + <)75( 
^ r) -(7« + <)7 5 (^) 



1 + / 
2 



7/?r 



2 



(10a) 
(10b) 



Possible contractions of VL a p are listed in the Appendix. From the Eqs. (3) and (8), Eq. (1) 



can be reexpressed in terms of and Q a/ j: 

V 4M A c i(l/2)^A 6 V3 L V 



(11) 
(12) 



A similar expression can be obtained for the spin-3/2 final states 
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(AMv',s')\T\A b (v : s)) 

y/AM Acl{3/2) M Ab 



t- T-i ( ( c ) i ( c ) ' ( c ) \ /3-p 

4^ a r - € e \Ti 'v a V/3 + T 2 V + ^ Q/3 J 7T 



(13) 



III. QCD SUM RULE EVALUATION 

As a starting point of QCD sum rule calculation, let us consider the interpolating field 
of heavy baryons. The heavy baryon current is generally expressed as 

3jA x ) = ^ T ^)CTj, P Tqi{x)}T\ P h k v {x) , (14) 

where i, j, k are the color indices, C is the charge conjugation matrix, and r is the isospin 
matrix while q(x) is a light quark field. Tjp and T' JP are some gamma matrices which 
describe the structure of the baryon with spin-parity J p . Usually Y and r' with least 
number of derivatives are used in the QCD sum rule method. The sum rules then have 
better convergence in the high energy region and often have better stability. For the ground 
state heavy baryon, we use ^1/2,+ — 7s, l"i/2 + = 1- ^ n ^ ne P rev i° us work [10], two kinds of 
interpolating fields are introduced to represent the excited heavy baryon. In this work, we 
find that only the interpolating field of transverse derivative is adequate for the analysis. 
Nonderivative interpolating field results in a vanishing perturbative contribution. The choice 
of T and T' with derivatives for the A^{ 2 and A 3 C { 2 is then 



Ti/2,- = (a + ty) 75 , r' 1/2j _ = ^75 , 



M 

r 3/2> _ = (a + ^) 75 , r 3/2; _ = ^(i^t + i%it) , (is) 

where a transverse vector is defined to be = N 1 — v^v ■ A, and M in Eq. (15) is some 
hadronic mass scale, a, b are arbitrary numbers between and 1. 
The baryonic decay constants in the HQET are defined as follows, 

(0bT/2,+ |A 6 ) = M v , (16a) 
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(16b) 
(16c) 



(ObT A -|Aci /2 ) = /i/2^ /2 , 

(o|^,_|Ai /2 > = ^h/^T" , 

where /i/ 2 and / 3 / 2 are equivalent since A^{ 2 and A^{ 2 belong to the same doublet with 
j[ = 1~. The QCD sum rule calculations give [9] 



f2 -A/T 
JA h e 



and [10] 

M / 1/2 e ' 



1 



20tt 4 



1 



du^e'^ + Uqq) 2 e 
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a 2 /8T 2 (ai s GG) 2 
+ 32tt 3 



(17) 



3iV c ! 



cu 7 (24a 2 + 406 2 )e 



+ ( a sGG) T ,^_ a 2 + ft 2) 
327T 3 



2tt 2 



4tt 2 



T M (3afc) . (18) 



{qq)T' 5 (lQab) - (qga ■ Gq}T' 3 ab - 

In the above equations, T^> are the Borel parameters and are the continuum thresholds, 
and N c = 3 is the color number. In the heavy quark limit, the mass parameters A and A' 
are defined as 



A' = M Aqi - m Q , A = M Aq - m Q . 



(19) 



The main point in QCD sum rules for the IW function is to study the analytic properties 
of the 3-point correlators, 

S hadron (u;, a/, y) ^ w V « n (cr) _ ^ ^ V + r 2 (c W + r s (c) <^) n<£> 



+e*(r< b) T;V + r 2 ( V^ + rf^)n^ 
Ei(u>,u>',y)=i*J ^ze^*"*^ 

'1 + f 



(20a) 



(20b) 



The variables k, k' denote residual "off-shell" momenta which are related to the momenta 
P of the heavy quark in the initial state and P' in the final state by k = P — rngv, 
k! = P' — m,Q/v', respectively. 
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The coefficient S(u;, u/, y)hadron in Eq. (20) is an analytic function in the "off-shell en- 
ergies" uj — v ■ k and uj' — v' • k' with discontinuities for positive values of these variables. 
It furthermore depends on the velocity transfer y = v ■ v', which is fixed at its physical 
region for the process under consideration. By saturating with physical intermediate states 
in HQET, one finds the hadronic representation of the correlators as following 

f 2 J* 

£hadron(^, v) = 7T, W7T V + hi g her resonances . (21) 

(A — uj J (A — uj) 

In obtaining the above expression the Dirac and Rartia-Schwinger spinor sums 

i + i 



(22) 



(23) 



8 = 1 

Af = *) = ( - nT + \iH) ^ , 

s=l ^ ' 

have been used, where gjf v = g^ v — v M t> v . 

In the quark-gluon language, uj', y) in Eq. (20) is written as 

poo n pclt (u v' v) 

E(u, u/,y)i3= / dvdv'-^- \; ,y) - + (subtraction) + S cond (cj, uj', y) 

2'2 J (u - u)(i/ - u') 

where the perturbative spectral density function p pert (z/, i/,y) and the condensate contribu- 
tion H cond are related to the calculation of the Feynman diagrams depicted in Fig. 1 . In Eq. 
(23), the 7-structures of spin-1/2 and 3/2 are the same as those in Eq. (20), respectively. 
Subleading IW functions, t- Q \ obtained from spin-1/2 and 3/2 are therefore identical. 
The six t-^ (Q = c,b, i = 1, 2, 3) are not independent. From the fact that 

id a (h®rh®) = hi c \tD a T + Tifi a )hf = (Av a - K'v'^Th^ , (24) 

Eq. (8) implies 

(r[ c) + r?)v a v p + (r 2 (c) + r 2 ( V^ + (r 3 (c) + r 3 % a/3 = (A.vp - A'^Ktfy) ■ (25) 

The above expression relates with rf^ as 

r[ c) + r? = M , (26a) 
r 2 (c) + r 2 (6) = -A'e , (26b) 
r 3 (c) + r 3 (6) = . (26c) 



Other relations are obtained from the equation of motion of the heavy quark, v ■ Dhv = 0: 



h$iv ■ *DrhW = {y T [ c) + r 2 (c) )rA„ = , (27a) 
h$nv ■ = r v ,T(r? + yr? + t^)a v = , (27b) 

From the above 5 equations in Eq. (26), (27), all the six subleading IW functions are reduced 
to only one independent form factor. We just pick up rf\y) = r(y), then others are 

r[ c) =M~r, (28a) 

4 C) = -yM + yr , (28b) 

r 3 (c) = y(yA - A')£ - (y 2 - l)r , (28c) 

r 2 (6) = (yA - A')£ - yr , (28d) 

rf = -y(yA - A')£ + (y 2 - l)r , (28e) 

Now that all the subleading IW functions are related to r(y), we have only to extract 
the coefficient of f^t^fi^p (or A^ a f Q ,f / 3r7 /3 for spin 3/2) in Eqs. (20) and (23). 

The QCD sum rule is obtained by equating the phenomenological and theoretical ex- 
pressions for S. In doing this the quark-hadron duality needs to be assumed to model the 
contributions of higher resonance part of Eq. (21). Generally speaking, the duality is to 
simulate the resonance contribution by the perturbative part above some thresholds u c and 
u' c , that is 



oo poo 



re,= / / ^ y > . (29) 



In the QCD sum rule analysis for B semileptonic decays into ground state D mesons, it 
was argued by Neubert in [20], and Blok and Shifman in [21] that the perturbative and the 
hadronic spectral densities can not be locally dual to each other, and therefore the necessary 
way to restore duality is to integrate the spectral densities over the "off-diagonal" variable 
v_ = \j^f\{y — vlS )l^i keeping the "diagonal" variable v + — {v + v')/2 fixed. It is in v + 
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that the quark-hadron duality is assumed for the integrated spectral densities. The same 
prescription shall be adopted in the following analysis. On the other hand, in order to 
suppress the contributions of higher resonance states a double Borel transformation in uj 
and u/ is performed to both sides of the sum rule, which introduces two Borel parameters 
T x and T 2 . 

Combining Eqs. (21), (23), our duality assumption and making the double Borel trans- 
formation, one obtains the sum rule for £(y) as follows; 



Mf 1/2 f: b e 



-A'/2T' -A/2T 



= 2 



y + l 

'jj' fyu) 
2T' D 2T 



V2 r wc(y) 



exp 



2T' 



2T 



_l_ R w z;cond 



(30) 



where i/ = i/+ + J^v-, u' = u+ - Jj+i^-, and 



i 



Cr = < 



(for spin 1/2) , 



(31) 



+e b (r 1 ( V^ + r 2 ( V^ + r 3 (6) ^jQ^ r) 

£v a r - e c (r[ c) v a vp + r 2 (c) w a ^ + r^g^^T 
+e b (rf ] v a vp + r 2 ( % a ^ + rf ] g a ^jT^ (for spin 3/2) . 

Now the remaining thing is to evaluate the relevant diagrams in Fig. I. The leading 
contributions are given in [12]. For the subleading corrections to the perturbative spectral 
density function p(u>,iv';y), we have 

p(u,,u';y) = B-'B-' 



a/3 



2v a v ,/3 f2cosh9A 3 B 3 e - e A 2 B 4 e e A 4 B 



sinh 6 



—1—0^)0(^)0(2^ - u 2 - u' 2 ) 
2 sinh 6 1 

4D2 



3!3! 



2!4! 



4!2! 



2 ^/3 / e 2e A 4 B 2 e -26 A 2 B A ^ 2A 3 B 3 
■ — —I — 77^ — 777777- 



sinh z #V 4!2! 2!4! 
from the perturbative diagram Fig. 1 (a), where 
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<X0 



A 3 B Z 
^!3T 



(32) 
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A = uj' - cue' 9 , B = ue 9 
e e =y+ ^y 2 -l . 



-LU' , 



(33a) 
(33b) 
(33c) 



For the condensate contributions we just give results when T' — T for simplicity; 



ibg a/3 fl 



a/3 



2tt 2 (1 + 2/) 2 
ibv a Q a p 



M{qq)T b -\(qgvGq)T\Ay + h/2) 



128 (qq)T 5 {3v + 2v'f 
+\(qg<y ■ Gq){(6y + 7/2)v? + (y - 3/2)^} 

2g a ^(2y 2 + 3y + 1) + (lOy + 6)v a v? + Ayv a v' p 
(y + l) 2 {y - 2){ - g a ? + 5v a (v + v'f) 



47T 2 (l + y) 3 

4 
3 

-' ' 12(1 • //) ;! ali 



-a 



a/3 



192tt 3 (1 + i/)' 

+24(y - l)v a v' p - 16(y + l)(y + 4)uV 

ia(a s GG)T 4 
~512vr 3 (l + |/) 4 



a/? 



-2(1 + y)# a/3 + 6t; c > + v 



(34a) 
(34b) 



(34c) 



Note that these results are from A^. If were the final state, VL a p would be replace by a 
proper 7-structure, leaving all the other things unchanged. 



IV. RESULTS AND DISCUSSIONS 

For the numerical analysis, the standard values of the condensates are used; 

(qq) = -(0.23 GeV) 3 , 
(aGG) = 0.04 GeV 4 , 
(qga ■ Gq) = m 2 (qq) , m\ = 0.8 GeV 2 . (35) 

There are many parameters engaged in the QCD sum rule calculations. The key point in the 
numerical analysis is to find a reasonable parameter space where the QCD sum rule results 
are stable. First, the continuum threshold u' c in (A') can differ from that in f\ b (A). 
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However, it is expected that the values of u c and uj' c would not be different significantly. 
This is because the mass difference A' — A is fairly small [10], A' — A ~ 0.2 GeV. Indeed 
the central values of them were close to each other in the sum rules analysis for /|(|) (A') 
and f\ b (A). One more thing to be noticed here is that the continuum threshold u c in 
Eq. (30) can be a function of y in general. But for simplicity, we take it to be a constant 
uj c {y) = w c = w[, = w in the numerical analysis. In this sense, we use only one constant 
continuum threshold throughout the analysis. An alternative choice of u c (y) = (1 +y)u>o/2y 
is suggested in Ref. [20]. We find that this choice yields almost no numerical differences. 
This is because the kinematically allowed region is very narrow around the zero recoil. 

Second, there are input parameters of a and b in the interpolating fields in Eq. (15). They 
are the parameters that generalize pseudoscalar or axial-vector nature of the light degrees 
of freedom (Ti/2,3/2 m Eq. 15). In Ref. [10], a particular choice of (a, b) = (1,0) gives the 
best stability for the mass parameter A'. We adopt the same choice of (a, b) = (1, 0) in the 
present analysis. 

Third, there are two Borel parameters 7\ and T 2 distinct in general, corresponding to uj 
and u/ in S(u;, a/, y), respectively. We have taken 7\ = T 2 in the analysis. In Ref. [16] for B 
into excited charmed meson transition, the authors found a 10% increase in the leading IW 
function at zero recoil when T 2 /Ti = 1.5 as compared to the case when T\ = T 2 . It seems 
quite reasonable for one to expect that in the case of heavy baryon, the numerical results 
should be similar for the small variations around T 2 /Ti = I. 

In short, we adopt the same parameters used in [10,12] where the mass parameter and 
the leading IW function are calculated. It makes sense because the observables involved are 
directly related to the subleading IW function r(y) through Eq. (30). 

In Fig. 2, r is plotted as a function of (y, T). Figure 3 shows the stability of r(y = 1) for 
the Borel parameter. The sum rule window is 

0.1 < T < 1.0 (GeV) . (36) 

The upper and lower bounds are fixed such that the pole contribution amounts to 50% while 
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the condensate one to 12%. One notes that the window given in Eq. (36) overlaps those 
obtained in the Refs. [9,10,12]. Of course, this reflects the self-consistency of the sum rule 
analysis. In Fig. 4, we present the shape of r(y) for a fixed Borel parameter. We found that 

r(y)=r(l)[l-p 2 (y-l)], 

r(l) = -1.277_g;£ GeV, for co = 1.4 ± 0.1 GeV , 
p 2 = 2.76^;^ , for uj = 1.4 ± 0.1 GeV . (37) 

V. SUMMARY 

Subleading contributions of Oil/mq) to the A b — > A cl weak form factors are important 
because some of the form factors do not survive at the heavy quark limit, and other remaining 
form factors vanish at zero recoil. Using the QCD sum rules, we calculate the subleading 
IW function r(y) which appears in the current matching in the HQET at (9(1/toq). We 
obtain r(y) given by 

r(y) = -1.27[1 - 2.7Q(y - 1)] GeV . (38) 

The best stability is attained when the continuum threshold u>o — 1-4 GeV. The parameter 
space for the analysis is the same as previous one for the leading IW function. The fact that 
by using the same set of parameters the present sum rule window for the mass parameter, 
leading and NLO IW function overlaps the previous ones ensures the self-consistenncy of 
the QCD sum rules. Our results can be applied directly to the decay mode A b — > A cl £u, 
along with the use of the previous LO IW function, but a complete analysis at 0(1/i71q) 
requires the information on another NLO contributions from the HQET Lagrangian. 
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APPENDIX A: CONTRACTIONS OF 



af3 



After a simple algebra, possible contractions for Q a/3 are given by 



v' a n a/3 = 


o, 




n + y 

V 2 


ap 


n + y 

V 2 


9 "a/3 — 


n + y 

V 2 




V 2 


u u a/3 


/i + y 

V 2 


» "a/3 


/i+y 

V 2 




V 2 




/i + y 

V 2 


= 


n + y 

V 2 




/i+y 
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FIGURE CAPTIONS 



Fig. 1 

Feynman diagrams for the three-point function with derivative interpolating fields. Double 
line denotes the heavy quark. 

Fig. 2 

Three dimensional plot of r as a function of y and T in units of GeV. The continuum 
threshold is chosen to be co c (y) = 1.4 GeV. 

Fig. 3 

t(1) as a function of the Borel parameter T. Each graph corresponds to lu = 
1.2, 1.3, 1.4, 1.5, 1.6 GeV, respectively, from the top. 

Fig. 4 

r(y) at a fixed Borel parameter T = 0.34. Each graph corresponds to uo = 
1.2, 1.3, 1.4, 1.5, 1.6 GeV, respectively, from the top. 
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FIGURES 




FIG. 2. 
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